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Active and reactive power in stochastic resonance for energy harvesting 

Madoka Kubotafl Ryo Takahashi, and Takashi Hikihara 
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A power allocation to active and reactive power in stochastic resonance is discussed for energy harvesting 
from mechanical noise. It is confirmed that active power can be increased at stochastic resonance, in the same 
way of the relationship between energy and phase at an appropriate setting in resonance. 


INTRODUCTION 


the following formula: 


Since noise appears everywhere in our surroundings, the 
energy conversion from noise to controlled motion is a key to 
come it to use. As an energy harvester from thermal noise, the 
molecular sized brownian ratchet was suggested, but this ma¬ 
chine was verified its impossibility Qll. On the other hand, 
stochastic resonance (SR) has been suggested as a method 
for energy harvesting from noise ffl. SR is a phenomena 
in which noise with moderate strength enhances SNR (Sig¬ 
nal Noise Ratio) id 71]. SR is possibly to apply in biological 
sensory orga ns M |9|], and the applications have expanded in 
medical use iflQL b rain image enhancement ifTlI] . and nano 
sized transistor (u 10. 

Here, we discuss SR for harvesting energy from white noise 
and the method to enhance the energy. This letter develops a 
concept of power factor correction as same as electrical sys¬ 
tems, and surveys a power allocation to active and reactive 
power. Finally it is clarified that active power of controlled 
motion is maximized at SR. 
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where R(f) is given as zero-mean white Gaussian noise of 
auto-correlation function: 


< R(f)R{t + At) >= 2ymkT6(At). 
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Where m denotes oscillator mass of a harvester, x displace¬ 
ment, my damping constant, h cos fit sinusoidal external 
force, U{x) bi-stable potential, k Bolztman constant, and T 
noise temperature. () implies time differential d/df and < > 
ensemble average operation. We take dissipative energy of 
controlled motion at frequency Q as harvested energy. 

Without loss of generality, we can introduce the previously 
obtained theoretical relationship from lEifiS. In particular, 
Q and cp are defined in 0141] and 1151] . 


SYSTEM AND POWER EQUATIONS 

Suggested energy harvesters by the method of SR possess 
bi-stable potential ffl. Here is assumed that dynamical 
equation of the SR harvesters are represented altogether by 



FIG. 1: Behaviors of amplitude Q and phase f as functions of kT, 
where m = 0.02, ym = 3.00, h = 0.20, and fl = 0.04. Both Q and (p 
show their single peaks. 
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Response of amplitude Q and phase tp for noise intensity 
kT are shown in Fig.[T] Qi is a relaxation rate defined by 
d^t//dx“|j:=±i ifTsIl . W/2 corresponds to Kramers rate 10 at 
Wo = ldx^\x=±\, tob - and potential barrier 

AU - 1/4. According to the reaction of Q in Fig.[Tl SR ap¬ 
pears around kT = 0.1. When SR appear^ the dissipative 
energy becomes maximum as confirmed in yHst]. The results 
explain that the external force frequency Q matches the time 
lag to overcome the potential barrier for the harvester. 

On the other hand, when a harmonic oscillator has phase lag 
7r/2 to sinusoidal external force, so called resonance appears, 
then dissipative energy becomes maximum. Similar phenom¬ 
ena is also confirmed in nonlinear oscillato rs ill71] . In Fig.[T] 
phase lag tp becomes maximum at SR lflol] illSll . However, it 
has not been accounted for the relationship between phase lag 
and dissipative energy. 
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ACTIVE AND REACTIVE POWER 


Here we focus on the relationship of phase lag (p and energy 
flow in system of Eqs. ([1]) and (|2]i. Since there are two external 
forces; sinusoidal force and noise, it is difficult to decide each 
contribution. For energy harvesting by SR, it is necessary to 
see the energy flow. From Eq.([T]), the following relationship 
is obtained. 

m d , , dt/ 

- < X > - -my < i > - < - X > 

2 df dx 

+ < X > h cos Qf+ < R(f)x > . (7) 


Here allocates each term to input, active and reactive power. 
Through the mechanical-electrical analogy IlSll . electrical 
voltage and current correspond to externally given mechani¬ 
cal force and the velocity of oscillator, respectively, so that the 
correspondence of electrical active power to mechanical dis¬ 
sipative power appears. At the same time, the reactive power 
can be explained as the mechanical energy flow. Generally 
active and reactive power are averaged over a period. On the 
other hand, instantaneous input, active, and reactive power are 
depicted as follows: 

Input power: 

< xh cos Qf > + < R(f)x > 

- -/ ^) h cos Qf + ykT, 
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Reactive power: 
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where the following two equations derived from the Fokker- 
Plank equation are substituted into Eq.(|7]i- 
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In Eqs.®, I®, and (fTOl i. those power are consisted of follow¬ 
ing Pi (i = 1,2,3): 


1 ldV\ 1 IdVV^ 

Pi =-(—) h cos Qf, P 2 = — {-^) 

ym \ dx / ym \ dx / 




Here we average Pi (i = 1,2,3) over a period 27r/Q, and ex¬ 
press them as Pi (i = 1,2,3). Fig.|2]shows them as functions 
of noise strength kT. P 3 is in proportion to kT. Pi and P 2 
reach thier peak around SR. 



FIG. 2: Pi (red solid line), P 2 (blue dotted line), P 3 (green broken 
line) as a function of kT calculated theoretically. 


Here we discuss the physical meaning of Pi (i = 1,2,3). 
The power sources of Pi and P 3 are obvious. Pi is input 
power from sinusoidai force h cos Qf, because < xh cos Qf >= 
Pi can be derived from Eq. (fT2l) . P 3 is input power from noise 
P(f), because < R{t)x >= P 3 can be derived from Eqs.® and 
@. P 3 corresponds to dissipative power consumed by local¬ 
ized small scale vibration. On the other hand, P 2 can not be 
derived directly from the given forces, and the power source 
is not clearly defined. As you can understand, P 2 is a term 
which cancels out by adding both active (Eq.®) and reactive 
powers (Eg. lfTOl) ). However, at first, P 2 includes frequency Q, 
since P 2 depends on potential shape V{x, f), as in Eqs.® and 
®. In addition, P 2 becomes maximum around SR as shown 
in Fig.|2l Therefore P 2 seems a power to vibrate over the po¬ 
tential barrier periodically. 

Next, energy flow is discussed based on the above power 
allocation. The input power P 3 from noise goes directly to ac¬ 
tive power and is consumed. Fig.|2]andEqs.®, ®, (fTOl i leads 
that, when SR appears. Pi contributes to the reactive power 
mostly. Then P 2 increases, since reactive power (= -P 2 + Pi) 
of Eg. lfTOl) has a limit. Consequently active power P 2 includ¬ 
ing Q also increases. And P 2 equals to the dissipative power 
of controlled motion. The special feature is that, when SR ap¬ 
pears, phase lag (p and also P 2 are maximized at an appropriate 
noise intensity. Hence, it is reasonable to adjust dissipative 
power by altering kT similar to power factor correction in a 
circuit. 


NUMERICAL ESTIMATION 


Figure[3]shows time dependence of displacement x at kP = 
0.03, 0.07, and 0.40. SR appears around kT - 0.07. Fig.|4] 
shows numerically calculated Pi,P 2 from Eqs.® and ®. 
Blue (gray) bullets at each kT are values of ten trials, and 
red (black) squares are ensenble averaged value from the ten 
trials. Pi and P 2 become maximum locally around kT - 0.07 
where SR appears. However the peaks are not clear as in 
Fig.El P 2 increases at kT - 0.5. This is because of a reduc¬ 
tion in calculation accuracy as we can see that the variance of 
blue (gray) bullets increases with noise intensity kT growth. 
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FIG. 3: Time changes of displacement x under the noise intensity 
kT = 0.03,0.07,0.40 with h and Q. kept constant. 
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FIG. 4: Numerically estimated Pi and P 2 under different noise 
strength. 


CONCLUSION 

In this letter, SR is investigated as a method for energy 
harvesting from noise. When SR appears, the phase lag 
and dissipative power of frequency Q are maximized at an 
appropriate noise intensity. The result suggests the possibility 


of maximizing dissipative energy of controlled motion as 
same as the analogous to power factor correction. 
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